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THE PLATEAU OF GAMMA-RAY BURST: HINT FOR THE SOLIDIFICATION OF QUARK MATTER? 

Shi Dai 1 , Lixin Li 2 , Renxin Xu 1 

ABSTRACT 

The origin of the shallow decay segment in gamma-ray burst's (GRB) early light curves remains a 
mystery, especially those cases with a long-lived plateau followed by an abrupt falloff. In this paper, 
we propose a mechanism to understand the origin of the abrupt falloff after plateau by considering 
solidification of newborn quark stars with latent heat released as energy injection to GRB afterglow. 
We estimate the total latent heat released during the phase transition of quark stars from liquid to 
solid states, to be order of ~ 10 51 ergs, which is comparable to the emission energy in the shallow 
decay segment. We also estimate the time scale of radiating the latent heat through thermal photon 
emission, and find that the time scale agrees with observations. Based on our estimation, we analyze 
the process of energy injection to GRB afterglow. We show that the steady latent heat of quark star 
phase transition would continuously inject into GRB afterglow in a form similar to that of a Poynting- 
flux-dominated outflow and naturally produce the shallow decay phase and the abrupt falloff after 
plateau. We conclude that the latent heat of quark star phase transition would be an important 
contribution to the shallow decay radiation in GRB afterglow, and would explain the general features 
of GRB light curves (including the plateau), if pulsar-like stars are really (solid) quark stars. 
Subject headings: 7-rays: bursts, X-rays, neutron stars, elementary particles 



1. INTRODUCTION 

NASA's broadband (gamma-ray, X-ray, UV & optical) 
Swift satellite's (Gehrels et al. 2004) successful launch 
and operation opens a brand new era in the observation 
of the gamma-ray burst (GRB) phenomenon and has rev- 
olutionized our understanding of GRBs in many aspects 
(Meszaros 2006; Zhang 2007). The prompt slewing ca- 
pability of Swift allows us to swiftly catch the very early 
signals following the GRB prompt emission and the pre- 
cise localizations make it possible for ground-based follow 
up observations of most bursts. A large number of well- 
sampled X-ray light curves from tens of seconds to days 
past the GRB triggers (Burrows et al. 2004; Liang et al. 
2007) have been accumulated through Swift's observa- 
tion, which provide us a good chance to investigate early 
afterglow and study GRB systematically. 

A canonical light curve of X-ray afterglows as revealed 
by Swift (Zhang et al. 2006; Nousek et al. 2006; OBrien 
et al. 2006) is composed of five parts: prompt gamma- 
ray phase with a tail, shallow decay phase, normal decay 
phase, jetlike decay phase, and erratic X-ray flares. The 
physical origins of these segments have been widely dis- 
cussed in the literature (Zhang et al. 2006; Nousek et 
al. 2006; Dai et al. 2006; Liang et al. 2007). In this pa- 
per, we mainly focus on the shallow decay phase and the 
subsequent phase. On the one hand, the physical origin 
of the shallow decay phase is still a mystery. Different 
models, such as the energy injection models (Zhang et 
al. 2006; Nousek et al. 2006; Panaitescu et al. 2006b), a 
combination of the GRB tail with delayed onset of the 
afterglow emission (Kobayashi & Zhang 2007), off-beam 
jets (Toma et al. 2006; Eichler & Granot 2006), precur- 
sor activity (Ioka et al. 2006), two-component jets (Gra- 
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not et al. 2006; Jin et al. 2007), two-component emission 
model (Yamazaki 2009), varying microphysical param- 
eters (Ioka et al. 2006; Panaitescu et al. 2006a; Fan & 
Piran 2006; Granot et al. 2006) and so on, are hard to 
differentiate among each other from the X-ray observa- 
tions (Zhang 2007). On the other hand, some puzzling 
facts related to the shallow decay phase are revealed by 
Swift's observations, which can not be explained by cur- 
rent models. Two interesting facts are that the break 
between the shallow and the normal decay segments 
in the X-ray light curve for some GRBs is chromatic 
(Panaitescu et al. 2006a; Fan & Piran 2006) and that the 
light curve of GRBs like GRB 070110 shows a long-lived 
plateau followed by an abrupt falloff (the decay slope is 
about —9, with time zero at the trigger) (Liang et al. 
2007). The former fact suggests that the optical and X- 
ray emission in the shallow decay phase may not be the 
same component, and the later fact indicates an internal 
origin of the X-ray plateaus and continuous operation of 
a long-term central engine. Systematically analysis of 
the Swift X-Ray Telescope data suggests that the phys- 
ical origin of the shallow decay phase is diverse (Liang 
et al. 2007) and we expect to get more information from 
it including the physics of dense material and possible 
gravitational wave (Alessandra & Meszaros 2009). 

We note that GRB central engine may relate to the 
physics of cold matter at supra-nuclear density, which is 
now one of the daunting challenges in particle physics. 
Cold quark matter is conjectured to be in a solid state 
at realistic baryon densities of compact stars (e.g., Xu 
2010, for a review), we are then considering latent heat 
of quark star phase transition from liquid to solid as en- 
ergy injection to GRB afterglow in order to understand 
the feature of plateau followed by an abrupt falloff in 
some GRB's light curve. Quark stars, as possible na- 
ture of pulsars, are likely to form in GRB, no matter in 
the process of high- mass star collapses or merge of binary 
neutron stars. In the solid quark star model, which is sue- 
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cessful to understand a variety of pulsars' observational 
features (Xu 2003), it is likely that as the temperature of 
quark star drops after GRB, a phase transition from liq- 
uid to solid would happen (Xu & Liang 2009). Since the 
quark star in phase transition would emit energy with 
constant temperature and solid quark star would cool 
very fast due to its low heat capacity (Yu & Xu 2009), 
the latent heat of this quark star phase transition not 
only provides a long-term steady central engine, but also 
naturally shows a abrupt cutoff when the phase tran- 
sition ends. In this paper, we show that not only the 
energy released during the phase transition and the time 
scale of radiating latent heat agree with observations, but 
the process of energy injection to GRB afterglow is also 
reasonable. Thus latent heat of quark star phase transi- 
tion as energy injection to GRB afterglow would be an 
prospective model to understand the physical origin of 
the shallow decay phase of GRB light curve, especially 
for the feature of abrupt falloff after plateau. 

In Section 2, we estimate the latent heat of quark star 
phase transition in solid quark star model. Section 3 
discusses the process of energy injection. We conclude 
the results in Section 4. 

2. LATENT HEAT OF QUARK STAR PHASE TRANSITION 

To estimate the latent heat of quark star phase tran- 
sition from liquid to solid, we need to know the state of 
cold quark matter and the interaction between quarks. 
However, due to the non-perturbative effect of the strong 
interaction between quarks at low energy scale and the 
many-body problem of vast assemblies of interacting par- 
ticles, we can not describe the state of cold quark matter 
from first principle up to now. Yet, it is phenomeno- 
logically conjectured that astrophysical cold quark mat- 
ter could be in a solid state, and a variety of observa- 
tional features, which may challenge us in the hadron 
star model, could be naturally understood in the solid 
quark star model (Xu 2003). Recent results of relativis- 
ts heavy ion collision experiments also show that the 
interaction between quarks is very strong in hot quark- 
gluon plasma (Shuryak 2009), then it is reasonable to 
conjecture that the interaction between quarks should 
be stronger in cold quark matter. The strong interaction 
may then make quarks grouped in clusters, and if the 
residual interaction between quark clusters is stronger 
than their kinetic energy, each quark cluster could be 
trapped in the potential well and cold quark matter will 
be in a solid state (Xu 2010). 

Considering that a single quark cluster inside a quark 
star is assumed to be colorless, just like each molecule in 
a bulk of inert gas is electric neutral, Lai and Xu in their 
recent work (Lai & Xu 2009) used Lennard- Jones poten- 
tial to describe the interaction between quark clusters in 
quark stars and got the equation of state of quark stars. 
The interaction is expressed as 

u(r) = Cf)% (1) 

where Uq is the depth of the potential and r can be 
considered as the range of interaction. 

Based on this solid quark star model, we can then esti- 
mate the latent heat of quark star phase transition from 
liquid to solid. Rather than performing difficult molec- 
ular dynamics simulations of crystallization, we would 



prefer estimating the latent heat in order of magnitude 
by analogy with inert gas and common substances since 
the quark clusters are non-relativistic and the interac- 
tion is similar to common substances. In Tablel below, 
we list the melting heat, heat of vaporization, potential 
and ratio of melting heat to potential of inert gas and 
some common substances. The melting heat corresponds 
to the latent heat, and for substances with known heat 
of sublimation, we equal the potential to heat of subli- 
mation, otherwise we equal the potential to the sum of 
melting heat and heat of vaporization. 

From the data, we can see that for most substances the 
ratio of melting heat to potential is between 0.1 to 0.01. 
Considering that the interaction between quark clusters 
is similar to inert gas and is relatively strong, we choose 
the ratio of potential to melting heat to be / « 0.1 <~ 
0.01 for estimation. Then based on the solid quark star 
model proposed by Lai, choosing Uq = lOOMeV (Lai & 
Xu 2009), we can estimate the energy released by each 
quark cluster in the liquid to solid phase transition as 

Cluster ~ fU » 1 ~ 10MeV. (2) 

For quark star of one solar mass, M Q w 2 x 10 33 g, the 
number of baryon is n = 10 57 , then the total energy 
released during the phase transition can be estimated as 

E = E clustcr n « 10 51 ~ 10 52 ergs. (3) 

This order of magnitude agrees with the typical energy 
released in the shallow decay phase of GRB, that is to 
say, the latent heat of quark star phase transition from 
liquid to solid is sufficient to produce the plateau. 

On the other hand, according to the Lindcmann law 
that a solid melts when the root-mean-square amplitude 
of atomic vibrations exceeds a certain fraction of the 
equilibrium nearest neighbor distance, we can estimate 
the temperature of quark star phase transition and fur- 
ther more the time scale of radiating latent heat. In 
Mohazzabi and Behroozi's work in 1987, they obtained 
the expression of the ratio of the root-mean-square am- 
plitude of atomic vibrations to the equilibrium nearest 
neighbor distance for inert gas, and found that the Lin- 
dcmann law being well consistent with experiments (Mo- 
hazzabi & Behroozi 1987). As for our estimation, we 
can consult the ratio of potential to heat, T = Uq/UT, 
for common substances, and then get the temperature 
of quark star phase transition by analogy. For one- 
component plasma, T rts 175 (DeWitt et al. 2001), for 
multi-component plasma, T « 233 (Horowitz et al. 2007). 
So choosing T w 200, U = lOOMeV, the temperature of 
quark star phase transition is around 0.5MeV. 

During the liquid to solid phase transition, the tem- 
perature of quark star would remain constant, and the 
latent heat would be released through thermal emission. 
Then the time scale of radiation can be estimated as 

t = aTMirR 2 ' (4) 

where, E = 10 51 ergs, a is the Stefan-Boltzman con- 
stant, R = 10km is the radius of quark star, and 
T ss 0.5MeV, we find that the time scale of radiation 
is t ~ 1000s, which agrees with observations of GRB 
afterglow plateau. 
After solidification, since the heat capacity of solid 
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quark stars is very low (Yu & Xu 2009), the central quark 
star would cool rapidly. The residual inner energy would 
be released almost immediately, thus an abrupt cutoff of 
energy injection to the afterglow would appear, and nat- 
urally result in an abrupt falloff after the shallow decay 
phase. In Figl, a schematic cooling curve of quark star 
is presented, which consists of three stages. The first 
stage could be an initial fast cooling stage due to the 
emission of neutrinos and photons at the very beginning 
of a quark star, when its initial temperature T could 
be much higher than lOMcV. When the temperature of 
quark star drops to the melting point T p , it could come 
to the second stage, the liquid to solid phase transition 
which last from t, to tf. At this stage, the temperature 
of quark star would remain constant and the latent heat 
would be released steadily which could provide a long- 
lived steady central engine. After phase transition, the 
born solid quark star would rapidly release its residual 
inner energy due to its low heat capacity, thus a steep 
falloff appears in the cooling curve corresponding to the 
abrupt falloff after the shallow decay segment. 

3. THE STANDARD FIREBALL REVISITED? 

In the standard GRB fireball model, it is assumed that 
a large amount of energy is instantly released through 
some explosion process. The energy drives some ma- 
terial to expand with a ultra-relativistic speed, which 
requires that the fireball outflow is low-baryon-loaded 
so that the total rest mass of the fireball is sufficiently 
small. The standard fireball is assumed to be highly non- 
uniform and in the extreme case is composed of many dis- 
tinct shells. Collision of different shells produces internal 
shock waves, which are thought to be responsible for the 
observed prompt emission. Collision of shells with the 
surrounding intermediate stellar material produces ex- 
ternal shock waves, which are thought to be responsible 
for the observed GRB afterglow emission. While for the 
afterglow flares, they are usually thought to be produced 
by later injection of energy through internal shock. 

After having estimated the energy released by quark 
star phase transition and the time scale of radiation, it 
is necessary for us to qualitatively discuss the process 
of energy injection to GRB afterglow. Generally speak- 
ing, the energy injection to afterglow could consist of 
some kinetic-energy (i.e., baryons) dominated shells or 
a Poynting-flux-dominated wind (Usov 1994; Meszaros 
& Rees 1997). In case of afterglow, we always con- 
sider an impulsive shell that is already heated during 
the shell-ISM interaction and that is collecting material 
from the ISM, and in the meantime also receives a large 
enough injection energy from a continuous Poynting-flux- 
dominated wind or a kinetic-energy dominated shell. As 
discussed in literatures (Zhang & Meszaros 2001, 2002), 
for Poynting-flux-dominated wind case in which pure en- 
ergy with negligible baryon loading is injected to fireball, 
no reverse shock is expected, and the injection signature 
is produced only by the forward shock emission. For 
kinetic-energy-dominated matter shells case, depending 
on whether the collision between the injected and the 
impulsive shell is mild or violent, the injection process is 
quite different. If the relative velocity between the collid- 
ing shells does not exceed a critical value defined by their 
energy ratio, the collision is mild, and the injection may 
be analogous to the Poynting-flux injection case. Other- 



wise, the injection is violent, and an additional pair of 
strong shocks will form at the discontinuity between two 
colliding shells which will greatly influence the injection 
signature. 

As for our energy injection model, the latent heat re- 
leased during quark star phase transition would produce 
a radiation-dominated fireball, whose luminosity is about 
L « 10 48 ergs/s according to our estimation. The optical 
depth can then be estimated as 

_ f p a T FD 2 

T77_ E?m e C 2 



•'^OOOMpc 7 v 10- 1:L crgs/cm 2n 10W ' v ; 

where, R = 10km is the radius of quark star. The optical 
depth is very large, so the energy flux from central quark 
star would form a pure radiation fireball consist of radi- 
ation and electron-positron pairs, which is similar to the 
fireball of prompt emission but continuous instead. We 
expect that this pure radiation fireball would expand and 
finally inject energy to the afterglow in form of Poynting 
flux. If we consider possible baryons in the space the 
fireball swept, we can estimate the velocity the matter 
shell could reach assuming that all the energy are con- 
verted to the kinetic energy of the baryons. Since in this 
late injection phase the baryon loading could be in prin- 
ciple much lower, we assume that the density of baryons 
are 10% of the original density in space. Then we can 
estimate as 

7 = E/Mc 2 . (6) 

For E w 10 51 ergs, we get 7 w 100. According to 
Zhang and Meszaros' work in 2002 (Zhang & Meszaros 
2002), for such kinetic-energy-dominated matter shells 
with 7 « 100 and L « 10 48 ergs/s, the collision between 
the injected and the impulsive shell is mild, and the in- 
jection may be analogous to the Poynting-flux injection 
case. 

4. CONCLUSIONS 

A possible physical origin of the shallow decay phase 
of GRB light curve is proposed in the solid quark star 
model. We suggest that quark stars may be born in 
GRB, and as the temperature of quark star drops, a 
phase transition from liquid to solid may occur in the 
quark star. The latent heat of phase transition would 
provide a long-lived steady energy injection to GRB af- 
terglow, since the temperature of central star would re- 
main constant during phase transition. When phase 
transition ends, an abrupt falloff after the plateau in the 
light curve would naturally appear which is hard to un- 
derstand by other central engine models. 

We estimate the latent heat of quark star phase tran- 
sition from liquid to solid based on the solid quark 
star model whose interaction between quark clusters is 
described by Lennard- Jones potential. The energy of 
~ 10 51 ergs and the radiation time scale of ~ 10 3 sec- 
onds agree with observations of the shallow decay phase. 
We also qualitatively discuss the process of energy injec- 
tion to afterglow, and show that the energy injection of 
phase transition would be in a form similar to that of 
the Poynting-flux-dominated outflow. Both the estima- 
tion and the injection process suggest that the idea of 
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considering latent heat of quark star phase transition as 
energy injection to GRB afterglow is rational. 
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Fig. 1. — A schematic cooling behavior of a new-born quark star. Three stages are shown: an initial cooling stage due to the emission 
of neutrinos and photons at the very beginning of a quark star with initial temperature To, a liquid to solid phase transition stage from 
time ti to tf with constant temperature T p , and, after solidification, a fast cooling stage because of solid quark star's low heat capacity. 
We focus on the duration from ti to tj in this paper. 
A 
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1MeV 




t 



TABLE 1 

The melting heat, heat of vaporization, potential and ratio of melting heat to potential of inert gas and some common 

substances (dean 1999). 





melting heat kcal/mol 


heat of vaporization kcal /mol 


potential kcal/mol 


melting heat /potential 


He 


0.0033 


0.0194 


2.2944 


0.0014 


Ne 


0.0801 


0.422 


9.1776 


0.0087 


Xe 


0.5495 


3.02 


12.6192 


0.0436 


Rn 


0.69 


4.01 


19.5024 


0.0354 












Al 


2.56 


69.5 


78 


0.0328 


Cs 


0.499 


16.198 


18.3 


0.0272 


Cu 


3.17 


72.74 


81 


0.0391 


Fc 


3.63 


83.68 


99.5 


0.0365 


Hg 


0.5486 


14.13 


14.65 


0.0375 


Na 


0.622 


23.285 


25.75 


0.0242 


Si 


12 


85.8 


107.7 


0.1114 


C 


25 




171.29 


0.1460 


CO 


0.2 


1.444 


1.644 


0.1217 


co 2 


1.99 




6.03 


0.3300 


H 2 


1.436 


9.717 


11.153 


0.1287 


H2O2 


2.987 


10.53 


12.34 


0.2421 


CaCl 2 


6.8 


56.2 


77.5 


0.0877 



